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The expected event rates for Zi dark matter for a variety of dark mat- 



■ r< ■ supergravity grand unified models. Radiative breaking constraints are im- 
p • plemented and effects of the heavy netural Higgs included as well as loop 
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^ ■ corrections to the neutral Higgs sector. The parameter space is restricted 
, so that the Zi relic density obeys 0.10 < f^^^/i^ < 0.35, consistent with the 
\-{ ' COBE data and astronomical determinations of the Hubble constant. It is 
found that the best detectors sensitive to coherrent Zi scattering (e.g. Pb) 
is about 5-10 more sensitive than those based on incoherrent spin depen- 
dent scattering (e.g. CaF). In general, the dark matter detectors are most 
sensistive to the large tan (3 and small nio and rrig sector of the parameter 
space. 
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1. INTRODUCTION 

There is much astronomical evidence that more than 90% of our Galaxy, 
and perhaps of the universe is made up of dark matter of unknown type. In 
galaxies, this matter has been detected by its gravitational effects on the 
motion of stars and gas clouds. A large number of candidates for dark mat- 
ter have been suggested both from astronomy and particle physics. In this 
paper we will limit our discussion to supersymmetry models with R parity, as 
they offer a natural candidate for dark matter, the lightest supersymmetric 
particle (LSP) which is absolutely stable. Thus the relic LSP left over from 
the big bang could be the dark matter present today. Further, in super- 
gravity GUT models, for almost all the parameter space of most models, the 
LSP is the lightest neutralino, the Z\. (The alternate possibility, that the 
sneutrino is lightest occurs only rarely.) Thus we will consider here only the 
Zi dark matter candidate, and do so within the framework of supergravity 
grand unification with radiative breaking. 

In this paper we discuss the expected event rates for a number of dark 
matter detectors using the following nuclei: "^i?e, ^°Ca ^^^2, ^^Ge -|-^^ Ge, 
^^Ga ^^^s, ^^Na^'^'^I and '^^'^Pb. The first two represent nuclei which are 
most sensitive to spin dependent incoherrent scattering of Zi by the nuclei, 
while the last four are increasingly sensitive to coherrent scattering. Pb could 
be a candidate for a superconducting detector. 

A great deal of work has already been done on the question of dark 
matter detector rates [1-7]. We present here an analysis over the entire 
SUSY parameter space with tan /? <20 which takes into account several 
important effects not generally treated before: 

• Radiative breaking. Almost all previous analysis has been done within 
the framework of the MSSM which does not include the constraints of 
radiative breaking of SU(2)xU(l). These constraints allow the deter- 
mination of //^ and niA [jJ- is the Hi — H2 Higgs mixing parameter, A 
is the CP odd Higgs boson) in terms of the other parameters. (Some 
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previous analyses have varied tjia arbitarily, obtaining spuriously large 
event rates.) 

• As pointed out in Refs. [6,7] the heavy neutral Higgs, H, can make an 
important contribution to the event rates. We have included this for the 
entire parameter space and find that the H contribution relative to the 
light Higgs, h, can range from 1/10 to 10 times as large. 

• As is well known, loop corrections to rrih are important due to the fact 
that the t quark is heavy [8]. We have also included the loop correction 
to a (the rotation angle arising in diagonalizing the h-H mass matrix). 
These actually cancel much of the effects of the loop corrections to m/j. 

• The COBE constraints on the Zi relic density are included. This strongly 
limits the region of SUSY parameter space that is allowed. In calculat- 
ing these relic density constraints it is essential to include the effects of 
the h and Z s-channel poles [9-11] for gluinos with mass 450 GeV. 
There are several effects we have not included here. Most noteworthy 

are that we have ommitted the possible WW, ZZ, Zh, hh final states in the 
Zi annihilation for the relic density calculation (which can occur when m^^ 
gets to the upper end of its allowed spectrum i.e. m^^^Mw and we have 
followed Refs. [12,13] in calculating the relic density. We estimate that this 
may lead to a (25-30)% error in the relic density, and since we have been 
reasonably generous in the allowed values for the relic density, we expect this 
will not significantly change our final conclusions. We also discuss below the 
sensitivity of the results to changes in the allowed region of Zi relic density. 

II. RELIC DENSITY CONSTRAINT 

The COBE data suggests that dark matter is a mix of cold dark matter, 
CDM, (which we are assuming here to be the relic Zi) and hot dark matter, 
HDM (possibly massive neutrinos) in the ratio of 2:1. In addition there may 
also be baryonic dark matter, B, (possibly brown dwarfs) of amount < 10% 
of the total. Defining ilj = pi/pc, where pi is the mass density of the i*'* 
constituent and pc — 3H'^/{8t:Gn) [H = Hubble constant, Gn = Newtonian 
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constant] is the critical mass density to close the universe, then the inflation- 
ary scenario requires Sfij = 1. A reasonable mix of matter is then — ^-^i 
^HDM — 0.3 and Qb — 0.1. What can be calculated theoretically is fi^^/i^ 
where h = H/ (100 km/s Mpc). Astronomical observations give h = 0.5-0.75. 
Thus we are lead to the estimate 

= 0.10-0.35 (1) 

Eq.(l) strongly resticts the allowed SUSY parameter space, and thus it 
is necessary to have a satisfactory method of calculating Q^^h'^. (We will 
discuss below the efi'ects of varying the maximum and minimum values of 
il^^/i^.) To do this, we use supergravity GUT models [14]. These models 
have the advantage of being consistent with the LEP results on unification 
of couplings at Mq ~ lO^^GeV [15], and generate naturally spontaneous 
breaking of supersymmetry in a hidden sector. In addition, by use of the 
renormalization group equations (RGE), the supersymmetry breaking inter- 
actions at Mq produce naturally spontaneous breaking of SU(2)xU(l) at the 
elctroweak scale Mz- In general, the low energy supersymmetry theory de- 
pends on only four parameters, ruo, rrig, At, tan/?, and the sign of fi. Here 
nio is the universal mass of all scalar fields at Mq, At is the t-quark cubic 
soft breaking parameter at the electroweak scale, and tan (3 — {H2) / {Hi) 
where (^^2,1) gives masses to the (up, down) quarks. 

The above may be contrasted with the MSSM (the formalism most dark 
matter calculations use) which possesses no theoretical mechanism for SUSY 
or SU(2)xU(l) breaking and is generally parameterized by 20 aribtrary con- 
stants. In the supergravity models, all properties of the 32 SUSY particles 
(masses, widths, cross sections, etc.) are determined in terms of the four ba- 
sic parameters and one sign. In particular, this means that tua and ji are so 
determined and are not free parameters (as usually assumed in the MSSM). 
Further, one finds throughout most of the parameter space the following 
(approximate) relations [16]: 

4 



while m/i<130GeF, m\ » m\ and tan (3 > 1. (Here, W\^2 are the two 
charginos and ^1,2,3,4 are the four neutrahnos). These relatfons will be im- 
portant in understanding the results below. 

The calculation of fi^^/i^ now proceeds in a standard manner. Using 
the RGE, we first express all SUSY masses and couplings in terms of the 
four basic parameters. This is done for the parameter space over the range 



ISOGeV <m~g< ITeV; lOOGeV <mo< ITeV; -2 < At/nio < 6; 2 < tan(3 < 20 

(3) 

with a mesh Aruo = 100 GeV, Arrig = 25 GeV,A{At/mo) = 0.5, and A{tanl3) 
= 2 or 4. We assume a top quark mass of mt = 167 GeV, and LEP and CDF 
bounds are imposed on the SUSY spectrum. The At range stated above 
exhauts the parameter space. Note that our analysis does not assume any- 
specific grand unification group but only that it is ai = (5/3)Q;y that unifies 
at Mg- in the early universe, the Zi is in equilibrium with quarks, leptons, 
etc. When the annihilation rate falls below the expansion rate, "freezeout" 
occurs at temperature Tf. The Zi then continues to annihilate via s-channel 
h and Z poles {Zi + Zi ^ h, Z ^ qq; ££; etc.) and t and u-channel squark 
and slepton poles. The relic density at present time is given by [13]: 



n,..-.4xio-(^)'(|l)'^ (4) 

where Nf is the effective number of degrees of freedom, (T^^/T^)^ is the 
reheating factor and 



J{^f) — I dx < av >; x — kT/rriz^ (5) 

Jo 

Here a is the annihilation cross section, v is the relative velocity and 
<> means thermal average. Since annihilation occurs non-relativistically. 



Xf 1/20, one may take the thermal average over a Boltzman distribution. 
However, as stressed in Refs. [9,10,11] one cannot generally make the non- 
relativisitic expansion av = a + bv'^ + ... due to the presence of the narrow h 
and Z s-channel poles. Thus calling flapprox the evaluation using the low v 
expansion, and 17 the rigorous result, we find for > that the relation 0.75 
< O,approx/^ < 1-25 is satisfied for only 35 % of the mesh points for rrig < 
450 GeV, but for almost 100 % for rrig > 450 GeV. The reason for this can 
be seen from Eq. (2). One is close to an s-channel pole when 2 m^^ ~ ^"^g 
is near ruh or Mz- Since m/i^l30 GeV, this cannot happen when mg^450 
GeV but one is usually somewhat near either an h or Z pole when rrig^ 450 
GeV. Thus a rigorous calculation is necessary for lower mass gluinos. 

The annihilation cross section a can be expressed in terms of the four 
basic parameters nio, rrig, At and tan /3. Using then Eq. (4) the region in 
parameter space obeying the COBE constraint of Eq. (1) can be determined. 

III. EVENT RATE CALCULATION 

Dark matter detectors see the incident Zi from effects of its scattering 
on quarks in the nuclei of the detector. This scattering proceeds through 
s-channel squark poles (Zi + q ^ q ^ Zi + q) and t-channel h, H and Z 
poles. These are some of the crossed diagrams to the annihilation diagrams 
appearing in the relic density analysis. Thus to a rough approximation, one 
may expect the event rate to be large when the annihilation cross section is 
large i.e. when fi^^/i^ is small. This makes results somewhat sensitive to 
where the lower bound on fi^^/i^ is set, and we will discuss this below. 

The scattering diagrams have been analysed by a number of people [1-7] , 
and we follow the analysis of Ref. [5].* One may represent the diagrams by 
the effective Lagrangian 

J^eff = {Xiri''xi)qr{A,PL + BqPR)q + {xiXi)Cgm,qq (6) 
* We include an extra factor of 4 in the cross section, due to the Majorana 
nature of the Zi , in agreement with Ref. [7] . 
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where q(x) is the quark field, ruq is its mass, and Xi{^) is the Zi field. Aq 
and Bq arise from the Z t-channel pole and q s-channel pole, and Cq from the 
h, H t-channel poles and q s-channel pole. Expressions are given for A,B,C 
in Ref. [5] . The first term of Eq. (6) give rise to spin dependent incoherrent 
scattering while the second term gives rise to coherrent scattering. There 
are several points to be made concerning the latter amplitude. In general, 
the Zi is a linear combination of two gauginos and two Higgsinos: 

Xi = aWs + PB+ 7#2 + SHi (7) 

The a, (3, 7, 6 can all be calculated in terms of the four basic parameters, and 
throughout the allowed part of the parameter space one finds 

P>a,S»-f (8) 
The coefficient Cq for the h and H poles is [17]: 
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where Fh = (a — ptan$w)i'ycosa + Ssina) and Fh = [a — |3tan9w){ls^na — 
5cosa). The tree value of a (the rotation angle that diagonalizes the h-H 
mass matrix) can be expressed in terms of the tree value of ruh [3] . Since loop 
corrections are large for the h particle, we have also included the loop correc- 
tions to a [8] in our calculation of Cq. Remarkably though, OLioop[{fnh)loop\ 
is generally quite close (within a few percent) to OLtree[{'^h)tree\- Thus a re- 
mains small i.e. a — 0{lQ~^rad). [Note, however, had one just inserted the 
loop correction to mu into the tree formula for a, one would have incorrectly 
obtained a large value for a, i.e. 5: = 0(1 rad.)!] One can now see why the 
H Higgs can make a significant contribution to Cq even though m^j » m^. 
For the d-quarks, the h term is reduced by a factor tan a relative to the H 
term. Further, the second fact in F^ is small, either because 7 is small or 
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sin a. is small. Thus for d-quarks, the H contribution can range from 1/10 to 
10 times the h contribution, depending on the point in the parameter space. 
[For u-quarks, the H term is generally small.] 
The total event rate is given by [5] 

R = {Rcoh+Rir.c)[PzJi^-^GeVcm-^)][< v^^ > / {320km/ s)][events/kgda] 

(10) 

where the coherrent and incoherrent rates are 

^'^oh = 16-/i— ^^-^210a I M,oh I' 

rrinr Mm „ „ 

^^^^ = ^^7W^ -^580AV(J+ l)C(r,p) | M,,e P (11) 

(Mat + m^J^ 

Here Mjv is the nuclear mass, C(^c/i)) C(^sp) are charge and spin form factor 
corrections, J is the nuclear spin and A is defined by < A?" | ES'j | N >= A < 
N \ J \ N > where Si is the spin of the z*'* nucleon. (A can be expressed 
in terms of the nuclear magnetic moment and nucleon g- factors.) Mcoh is 
proportional to Cq and Mine is proportional to Aq — Bq, explicit formulae 
being given in Ref. [5]. 

IV. RESULTS 

Eq. (11) allows one to divide dark matter detectors into two categories: 
those sensitive to the incoherrent (spin dependent) scattering due to a large 
value of A^J(J+1), and those sensitive to the coherrent scattering. Examples 
of "incohererent detectors" are ^He and ^°Ca ^^Fa with CaFa the most 
sensitive detector. Eqs. (11) show that Rcoh ~ Mjy and Rinc ~ '^/Mn 
for large Mn, the additional factor in Rcoh arising from the niq factor 
in Eq. (6), i.e. roughly speaking, all the quarks add coherrently to yield 
a Mn factor in the amplitude. The remaining detectors considered here, 
^^Ge Ge, '^^Ga'^^As, ^^Na ^^^j and ^'^'^Pb are all of the "coherrent" 
type with Pb being the most sensitive since it is heaviest. 
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The dependence of the expected event rate on the supergravity GUT 
parameters is fairly comphcated as each parameter enters in several places 
and the constraint Eq. (1) on fi^^/i^ limits the parameter space. One can, 
however, get a qualitative picture of the parameter dependence by studying 
several characteristic examples . Fig. 1 shows that R decreases rapidly with 
rrig (mainly because the Zi becomes more Bino-like). It also shows that R 
is larger for larger tan (3. (See e.g. the 1 /cos /? factor in the denominator of 
the d-quark part of Eq. (9); the 1/sin (3 factor for the u-quark part never 
gets exceptionally large since tan /3 > 1 in the radiative breaking scenario). 
Finally we note that R[Pb] is 5-10 times larger than R[CaF2] which is also a 
general feature. The tan /? dependence is shown more explicitly in Fig. 2 for 
the Nal and Ge detectors. (The three examples were chosen so the 0, h? is 
roughly the same at each tan j3 along each graph). The Nal curve lies higher 
than the Ge one for each pair since I is heavier than ^^Ge. 

In general, the event rate drops with increasing ruo, as one would expect 
since the squark mass increases with mo, reducing the eflFect of the s-channel 
squark pole. (There are additional effects, however, as nio also enters in the 
radiative breaking equations, effecting the size of //.) Fig. 3 illustrates the 
general behavior for several of the detectors. The coherrent detectors, Pb, 
Nal, Ge, scale almost exactly by their atomic numbers. (Fig. 3 also exhibits 
one of the few regions of parameter space where the CaF2 detector lies above 
the Pb detector.) 

Fig. 4 exhibits the maximum and minimum event rates for the Pb de- 
tector (the most sensitive of the coherrent detectors) and the CaF2 detector 
(the most sensitive of the incoherrent detectors) as a function of At, as all 
other parameters are varied over the entire space. One sees that generally a 
Pb detector will be a factor of 5-10 times more sensitive than a CaF2 detec- 
tor. Other coherrent detectors have event rates that scale with the Pb curve 
in proportion to their atomic number while the ^He has event rates a factor 
of 3 smaller than CaF. 
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The above analysis has been done with fi^^/i^ obeying the bounds of 
Eq. (1). We discuss now the eflFect of varying these upper and lower limits. 
As mentioned in Sec. Ill, the event rate R rises with decreasing fi^^/i^ 
and this rise is rapid near Oj^^h^ — 0.1. Further, the maximum value of R 
occurs when rrig is near its minimum value. However, by the scaling relations 
Eq. (2) this can force < A5GeV, and hence such parameter points are 
excluded by the LEP bounds. This is what causes the sharp peaks in Fig. 
4, which occur when lies just above the LEP cut. If, for example, one 

increases the lower bounds on fi^^/i^ to 0.15, one finds that the maximum 
event rates follow the curves of Fig. 4 with the peaks cut off. 

The upper bound on fi^^/i^ determines the minimum event rates. This 
is because the minimum rates occur when rrig takes on its largest value. As 
nig increases, so docs m^^ by Eq. (2). The Zi annihilation cross section then 
drops and ^Zi^'^ rises. The upper bound of Eq. (1) on fi^^/i^ is found to 
occur when {rrig) Max — 750 GeV. If one reduces the upper bound on fi^^/i^ to 
0.2 (which is consistent with the inflationary scenario which prefers h ~ 0.5). 
Then the maximum value of rrig is reduced to * {mg)Max — 400 GeV. This 
then increases the minimum event rate curves of Fig. 4 by about a factor of 
10. 

V. DETECTION POSSIBILITIES 

The above discussion has analysed the expected event rates for a variety 
of dark matter detectors over the range of parameters of supergravity GUT 
models. These detectors are most sensistive to the region of parameter space 
where tan (3 is large and rrio and rrig are small. Two types of detectors were 
noted: those with nuclei most sensitive to the spin dependent incoherrent 
scattering of the Zi (e.g. CaF2), and those most sensitive to coherrent 
scattering (e.g. Pb). In general, the best of the coherrent scatters are more 
sensitive than the incoherrent scatterers by a factor of 5-10. The coherrent 

* Such a low mass gluino could make it accessible to detection at the 
Tevatron. 
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scatterer event rates scale approximately with their atomic number. 

Dark matter detectors currently being built plan to obtain a sensitivity 
to signals with R^O.l events/kg da. Future detectors may be able to obtain 
a sensitivity of -R^O.Ol events/kg da. From Fig. 4 one sees that the present 
detectors will be able to study a small fraction of the total event rate, partic- 
ularly the large tan /? region. A more sizeable portion could become available 
to the next generation of detectors. However, it would appear that there will 
be regions of parameter space inaccessible to these types of dark matter de- 
tectors. In spite of this, the experimental study of even only a small part 
of the parameter space is of real importance, as such results, combined with 
other experiments (e.g. the 6 — > s-l-7 decay), can together significantly limit 
the allowed parameter space of supergravity grand unification models. Thus 
recently it has been demonstrated that the experimental limits on 6 — * s 4- 7 
from CLEO do indeed affect relic density analyses [18]. An analysis of the 
event rates with the inclusion of the CLEO constraint will be given elsewhere 
[19]. 
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FIGURE CAPTIONS 

Fig. 1 Event rates R vs nig for Pb (solid curve and CaF (dashed curve) de- 
tectors, with At /mo = 2.0, nio = 100 GeV and // <0. The upper hne 
of each curve is for tan (3 = 20, the lower line for tan /3 = 6. The gap 
between the short branches of the curves (at rrig = 225 GeV) and the 
main branches are regions where ^Zi^"^ ^ 

Fig. 2 Event rates R vs tan /3 for Nal and GeV detectors for nig = 275 GeV, 
fj, > 0. The solid curve is for At/nio = 0.0, nio = 200 GeV, the dashed 
curve for At /nio =0.5, nio = 300 GeV, the dash-dot curve for At /nio = 
1.0, nio = 200 GeV. The upper curve of each pair of curves is for Nal, 
the lower for Ge. 

Fig. 3 Event rates R vs ruo for Af/nio = 0.5, tan /3 = 8, nig= 300 GeV, > 
0. The dashed curve is for CaF2, and solid curves in descending order 
for Pb, Nal and Ge. 

Fig. 4 Maximum and minimum event rates vs. At /nio for Pb (solid) and CaF 
(dashed) detectors, for /j, < 0. All other parameters are varied over the 
whole space to obtain the maximum and minimum event rates. 
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